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Abstract 

In models with an enhanced coupling of the Higgs boson to the bottom quark, the dominant 
production mechanism in hadronic collisions is often the partonic sub-process, bg — > bH. We 
derive the weak corrections to this process and show that they can be accurately approximated 
by an "Improved Born Approximation". At the Tevatron, these corrections are negligible and are 
dwarfed by PDF and scale uncertainties for Mh < 200 GeV. At the LHC, the weak corrections 
are small for Mh < 500 GeV. For large Higgs boson masses, the corrections become significant 
and are ~ 18% for M H ~ 1 TeV at ^ = 10 TeV. 
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I. INTRODUCTION 



The search for the Higgs boson is one of the most important tasks for both the Fermilab 
Tevatron and the CERN Large Hadron Collider. The Standard Model requires a single 
scalar Higgs boson, with well defined properties except for its mass. In this case, the Higgs 
boson will be discovered at either the Tevatron or the LHC, with the discovery channel 
depending strongly on the Higgs boson mass. In the Standard Model, the production of 
a Higgs boson in association with b quarks is never a discovery channel due to the small 
b quark- Higgs boson Yukawa coupling. However, in non-standard models of electroweak 
symmetry breaking with a light Higgs boson, the coupling of the Higgs boson to the b quark 
is often enhanced and the channels bb — > H and gb — > bH become important modes 1-9]. 
A familiar example of such a model is the MSSM with large tan where enhancements by 
orders of mag nitude over the Standard Model prediction are possible in some parameter 



regions 



10|-ll3|. 



The hadronic production rate for the associated production of a Higgs boson and a b 
quark is well understood lHSj. The calculation can be done in either a 4- flavor or a 5- flavor 
number parton distribution function (PDF) scheme, which represent different orderings of 
perturbation theory. In the 4- flavor number PDF scheme, the lowest order processes for 
producing a Higgs boson and a b quark are gg — > bbH and qq — > bbH. Alternatively, in the 



arge logarithms of the form 
. In this scheme, the 



5- flavor number scheme, the b quark is treated as a parton and 
ln(^-) are absorbed into b quark parton distribution functions (14 
lowest order process for producing a Higgs boson in association with b quarks is bb — > H 
when no b quarks are tagged in the final state, and bg — > bH when a single outgoing b quark 
is tagged. Within the uncertainties, the 4— and 5— flavor number schemes give equivalent 
results for the NLO QCD corrected rate for associated b quark-Higgs production^. 

We work in the 5— flavor number scheme for simplicity and consider the associated 
production of a b quark and a Higgs boson. The bb — >• H rate is known to NNLO QCD 16], 
along with the full electroweak and SUSY QCD corrections 13, llSj. When an outgoing b 
quark is tagged, the rate is lower, but the background is significantly reduced, making this an 
important channel. Both the NLO QCD [l- 6 [and the SUSY QCD (SQCD) corrections from 

191 ] are known for bH production. Furthermore, 



gluino-squark loops in the case of the MSSM 



the Tevatron experiments have produced limits on bH production in the MSSM[2Cj,|2l|] which 
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can be interpreted in terms of the fundamental properties of the model. 

In this paper, we compute the Standard Model weak corrections to the bg — > bH process 
and compare them with the scale and PDF uncertainties of the NLO QCD corrected rates. 
We also compare our results with an approximation where the dominant corrections arise 
from the on-shell corrections to the bbH vertex (Improved Born Approximation). Section II 
contains the theoretical framework for the weak corrections. We retain the effects of a non- 
zero b quark mass everywhere. Numerical results are given in Section III and conclusions in 
Section IV. 



II. THEORETICAL FRAMEWORK 



In this paper, we consider the Standard Model process of associated b quark- Higgs boson 
production. Our results can be generalized in a straightforward manner to models with non- 
standard b quark - Higgs boson couplings. The tree level coupling of a b quark to a Standard 
Model Higgs boson, H, is given by 

Lyuk = —QboboboHo , (1) 

where the subscript, '0', denotes the unrenormalized quantity and 

mbo e f . 

im Wq s Wo 

We work in an on-shell scheme where the weak mixing angle is a derived quantity and is 
defined in terms of the physical gauge boson masses, 

sin 2 6 w = s 2 w = l-^. (3) 

z 

The lowest order Feynman diagrams for the process b(pi) + g(qi) — > bfa) + H(p 2 ) are 
shown in Fig. HJ The resulting Born cross section isjo], 



dajbg ^ bH)(fj, R ) 1 a s (m) - ( , 2 / M 4 H + u 2 2m 2 

dt {s-mlY 24 9b[flR) \ Sl tt + s\t\ L 



4u*isi + M 2 H (M 2 H -%i" 2 



' 2 f 2 
h 1 l 1 



M 2 H -u , (4) 



where si = (pi + q\) 2 — m 2 = s — mf, t\ = (pi — P2) 2 — ml, and u = (pi — q2) 2 , are the 
usual Mandelstam variables and the scale hr is the arbitrary renormalization scale. (The 
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FIG. 1: Lowest order Feynman diagrams for the process bg — > bH. 



cross section for the charge conjugate process, bg — > bH, is identical to Eq. HI) In the limit 
nib — > 0, the tree level contribution to bg — > bH vanishes and the first non-zero contributions 



are a subset of the 1- 
considered in Refs. 



22 



oop 



amplitudes computed in this work. The rrib — > limit has been 



23] and we will comment on the numerical effects of this limit in 
Section III. In this paper, however, we keep mj nonzero everywhere. 

In Eq. HI the Yukawa coupling, (^(^wr), is expressed in terms of the 1-loop renormalization 
group improved running MS mass for the b quark, Wib(iip>). For the decay H — > bb, the 



0(a s log(^-)) contributions can be absorbed into TrTh(Mu) [24i. [25[, motivating our use of 
the running mass. The O(as) NLO predictions for the production process, bg — > bH, 
however, depend sensitively on this choice for <7&|lj]. 

A. Renormalization 

As input parameters in the electroweak sector, we take a(0), Mz, and G^, along with 
the Higgs boson and fermion masses. The W mass is then a derived quantity. At tree level, 

W0 V2G, (1 - M^ /Mf ) • 1 1 

The gauge boson 2 -point functions are defined as, 

^xy(p 2 ) = 9^xy(p 2 ) + P^Bxvip 2 ) , (6) 

where XY = WW, Z 77 and jZ. Analytic results for the Standard Model gauge boson 
and Higgs contributions can be found in Refs. j^f], l^ 1 and for the fermion contributions 



1 A convenient compilation of the gauge boson and Higgs contributions to the gauge boson 2-point functions 
employing our conventions can be found in the appendix of Ref. (2^ . 
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in the appendix of Ref. 29(. The gauge boson mass renormalizations (My = Mw, Mz) are 



defined by the on-shell condition, 



M-M Hl + 5)^( 1 + «). m 



The electromagnetic charge renormalization is determined from Thompson scattering by 2 , 

/ 5a \ 

a = a 1 H 

V a J 

where the contribution from gauge bosons, leptons, and the top quark is given by IT 77 (0) = 
(9II 77 (p 2 ) / 'dp 2 | p 2 =0 . The light fermions contribute large logarithms to II' (0) which we 



77 

estimate by II' (0) \ had = iT 77 (Mf )/Mf + Aa 5 \ had , with Aa 5 \ had = .02761 ± .0036|30|. 



The Fermi constant is found from muon decay and, by definition, QED corrections 



are included in the measured value 



31 



32|, 



where, 



5G^ U ww (ti) 



a ( 7-As 2 Wl . 2 A 2 II 7Z /in , 

5 y _ B = --—2- 6 + 2 w log 4 + . 10 

4vrs^ V 2s 2 w J s w c w M| 

In the second line of Eq. [101 Sv-b is the contribution from vertex and box diagrams. For 
consistency, the W mass in this equation should be found from the tree level prediction, Eq. 

El 

The weak mixing angle renormalization is defined, 

s 2 wo = s 2 w (l + S 4 L )- (11) 



'it- 



Using the on-shell definition of Eq. El 



5s"w Cw 

( II W (M&) , Tlzz{M 2 ) , 

s 2 s 2 I M 2 M 2 ' ' 

6 iy *w V lvl w lvl z 



This relation effectively defines our convention for the sign of Ii 7 z. 
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The fermion renormalization proceeds in the usual manner 3 , 



m b0 = m b + 5m b 



°L,R 



5Z 



L,R 



>L,R ■ 



(13) 



The one-loop b quark self- energy is, 



Y, b {k 2 ) = k 



PA ^(k 2 ) + Re(5Z L ) )+Pr( Z«(k 2 ) + Re(5Z R 



-m b 



5mb v n.2\ i 5z l + 8Z ] r 5Z R + SZ[ 

-^;- Es{k) + — 2 — Pl + — 2 — Pr 



(14) 



Imposing the on-shell conditions, 
5m b 



m b 



Re{ EgOmi) + 



S^K) + S«(m 



V\"%j 



5Z L = -RetX^mt) -m 2 b [Zt(m 2 b )' + Z*(m 2 b )' + 2Z s (m, 



,2\/ 



8Z R = -Re{Y>*{m 2 b ) - m 2 b (^{m 2 b )' + S«(m fe 2 )' + 2E s (m^' 



(15) 



where S(m 2 ) / = dH(p 2 )/dp 2 \ p 2 =m 2. Analytic results for the expressions in Eq. [13] are given 
in Refs. flH- 



The Yukawa coupling renormalization is defined as, 



9b = 9b[l + 



(16) 



Combining the above results, the one-loop electroweak counterterm corresponding to Eq. [T] 
is then, 



Lyuk = ~9b (l + for J bbH 

5g b 5Z H 5Z L + 5Z R 

Oct — 

9b 2 2 



(17) 



where, 



$9b 
9b 



5m b 1 5G 



i-i 



m b 2 G M 
5m b 5Mw Se 
m b Mw e 



5s 



w 



s w 



(18) 



b L , R = (1 =F 75 )6 and P L , fl = (1 T 7s)/2. 
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FIG. 2: (a) Virtual QED corrections to the decay H — > bb, (b) Real photon emission contributions 
to H — > blyy, and (c) Feynman diagram contributing to the QED counterterms for H — > bb. 



Finally, we need the 1-loop corrected value for the W mass [3 31] . 



2 7ra 1 

W = V2G^ W 1 - Ar 




where, 



Ar = IT' ( )-2^^P + 



c w n 7Z (o) iw(o) - u ww (m^) 



ln 1 + Ar)), (19) 



77V 7 «w Mf M 2 , 



,2 



c^iWMf) imm,|.)\ W 6 + I^ log((4) , j 



s 2 V M 2 M 2 4tt<? 2 V 2s 2 

B. Electroweak corrections to H — > 66 

The calculation of the electroweak corrections to gb — > bH has many common features 
with that for H ^ bb and so we review the decay process briefly. This discussion will also 
serve to make clear our separation of the QED and weak contributions. The electroweak cor- 
rections contain both pure QED photonic contributions and the remaining weak corrections 



and the two contributions are separately gauge invariant 
The corrections to H — > bb can be parameterized as, 



34) 



T(H 66) = T(H bb) (l + A b ™ D + Agf D + Af&) , (21) 



r(H -> bb) 
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0.15 ■ 



I 1 1 1 I 

bbH 



'QED 
bbH 

'WK 



0.1 ■ 



0.05 



I 

I 



-0.05 1 



2(10 



I i i i I i i i I i i i I i i i I i i i I 
300 400 500 600 700 800 
M H (GeV) 



900 1000 



FIG. 3: Pure QED (black solid) and weak (red dashed) corrections to the decay H — )• bb as defined 
in Eq. EU 



where T(H — > bb) Q is the tree level result, but evaluated with a running b quark mass as 
described above. The QCD corrections are known to 0(a 2 ), including all top quark mass 



effects 
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36]. 



The QED corrections are found from the virtual photon diagram of Fig [2^, the real 
photon emission diagrams of Fig. [2b, and the counterterms derived from Fig. [2b, 



a bbH _ \bbH , a bbH , a bbH 

^QED — L ^Virt,QED ' L ^CT,QED "r L ^Real,QED 



(22) 



The virtual photon contribution in N = 4 — 2e dimensions is, 



a bbH 

L *Virt,QED 



a 



Q] 



71 



Li? 



1-/3 



2\ e 



+ 



Lin 



1 + 2 

P 

1 + P 



- X -L 

2e 



L 2 7T 2 



4 3 
1-0 2 

P 



+ 

L + 3 



a) ""'v 2/3 y j v p r ' ~j ' (23) 

where = yl — 4m 2 / M H and L = log((l — p)/(l + 0)). This is in agreement with Refs. 
371 . l38| with the replacement |a s — > Q\a and Qb= —\- 

To find the counterterms, we need the photonic contributions to Eqs. H7] and [18j By 
definition, <5G M contains only the weak contributions. Similarly, the Higgs boson self-energy 
does not receive contributions from diagrams with photons in the internal loops. Thus, we 



have to separate the QED contributions to Srrib, 5Zl and 5Zr which come only from Fig. 
[2b . The results are well known, 



sz L 

,QED — 

5Z R 

,QED 

f 5m b \ 

V m b J QED 

The QED counterterm is then, 

AcT,QED = 



ft 



4vr 6 V 



mi 



~—Q b 
-Qi\~ - 



An 



m 



r(i + e 
rfi + e) 



+ 4 



+ 4 



Sm b 5Z L 



5Z 



R 



--QI 

7T 



{ Ann* 



r(i 



Finally, we need the real photon emission contributions from Fig [2b, 



*Real,QED 



7T V m fe / 



1 + 



1 + 

2/3 



L 



+ F{M H ,(3) 



37- 



(24) 



(25) 



(26) 



where F(Mh-/3) is finite and an analytic form can be found in Refs. |34J. 

The QED contributions enumerated above are recognized by the explicit over-all factors 
of Of. The remaining weak corrections to H — ?■ bb, A^/^, are given analytically in Refs. 



34 



and are found from diagrams with W's, Z's, and Goldstone bosons, along with top 
quark contributions. In Fig. [3] we plot the QED and weak contributions to T(H — > bb). The 
QED corrections are always O(10~ 3 ) and can safely be neglected for all practical purposes 4 . 



C. One-Loop Corrections 

The one loop weak corrections to the process bg — > bH consist of self energy, vertex, 
and box diagrams (Figs. HH6]), along with the counterterms given explicitly in Eq. [T71 Our 
results can be expressed as, 

a(bg bH) NLO = a(bg bH) (l + A QCD + A QED + A WK ^j , (27) 

where o"o is the Born cross section of Eq. H] evaluated with the 1— loop renormalization 
group improved value for 7^. 

The purely photonic QED corrections consist of vertex and box contributions from inter- 
nal photons along with the corresponding counterterms, real radiation from bg — > bH'y, and 



The spikes at the W + W and ZZ thresholds are softened if the complex mass scheme is employed 4Cll41|. 
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FIG. 4: Self energy diagrams contributing to the weak corrections to bg — > bH. 



the process involving photons in the initial state, 76 — > bgH. The 0(a) QED corrections to 
bg — > bH can be found from the corresponding QCD corrections by making the substitution 
|a s — > aQl jl- 6 1. However, evaluating the process 76 — > bgH requires the use of a PDF set 
which includes initial state photons. 5 This contribution is expected to be quite small since 
potentially large logarithms from initial state collinear photon emission are absorbed into 



the PDFs. We further note that the Q 



the corresponding decay H — > bb 34 



]D contributions to the bb — )■ H process and to 
discussed above, are known to be less than 1 %. As 
this is considerably smaller than the PDF and scale uncertainties which we present in the 
next section, we do not provide numerical results for the pure QED corrections to bg — > bH, 
but evaluate only the weak corrections. 

The Feynman diagrams are generated using FeynArts 



43J and the interference with 



the tree level amplitude is evaluated numerically in Feynman gauge using FormCalc and 
Loop Tools 44| . We retain a non-zero bottom quark mass everywhere. 



D. Large Higgs Mass Limit 



The contributions to the weak corrections in the large Higgs mass limit can be easily found 
and provide a check of our results. The large Higgs mass limit for the process bg — > bH is 
obtained by noting that the triangle and box diagrams shown in Figs. H] -|6] are of O 



The most modern set of PDFs which include initial state photons are the MRST2004qed PDFs [42 1. 
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FIG. 5: Vertex diagrams contributing to the weak corrections to bg — > bH. 



relative to the tree level amplitude. The only contributions which are enhanced by factors 
of — come from the renormalization of the bbH vertex, Eq. [TJ In the large Mh limit, 

Lyuk -> 7777- ( z ^ , 2 \bbH + terms not enhanced by ^f- . (28) 
2M W V 1 + SM^ /Mfa J v 2 



11 




FIG. 6: Box diagrams contributing to the weak corrections to bg — > bH. 



The large M# limit of the Higgs wavefunction renormalization is 45|, |46[, 

Z H = 1 



1 Ml, r 
4- 6-ttV3 



167T 2 V 2 

Similarly, the W mass renormalization receives a contribution proportional to Mg, 

5M 2 , 1 ME 



32tT 2 V 2 ' 

The leading M# corrections to bH production are therefore 



45|, 



a(bg^bH) 



H 



■a(6p -> Mf) t 



1 + 



1 M 



li 



32tt 2 w 2 



13 - 2nV3 



a(bg -> MT)o 



1 + A™^°° )a(bg^bH) . 
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FIG. 7: Lowest order and NLO QCD results for pp — > b(b)HX at the Tevatron with y/s = 1.96 TeV 
with yfs = 1.96 GeV, ]P T > 20 GeV, | t/j, |< 2, and AR > .4. The renormalization/factorization 
scales are set equal to \i. 

III. RESULTS 

A. Numerical Results 

For our numerical studies, we use the following inputs: 

a = 1/137.03599911 
G M = 1.16637 x 1(T 5 GeV~ 2 
M z = 91.1875 GeV 

M t = 173.1 GeV . (32) 

We set the CKM mixing matrix to unity. For the pole mass of the b quark, we take Mb = 
4.25 GeV. We use CTEQ6.6 PDFsj^ and vary the renormalization/factorization scales 
f 1 — — \i-F from M H /2 to 2M H in the total cross section results. 
Our results are expressed as, 

a(bg ->• bH) NLO (fi) = a(bg ->■ bH) (fi) (l + A QC d(v) + ^qed(^) + &wk(^ , (33) 

13 
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FIG. 8: PDF uncertainties for pp -> b(b)HX at the Tevatron with y/s = 1.96 TeV p b T > 20 GeV, 
| rj h |< 2, AR > .4, and /x = M H /2. The solid line is a NLO (CTEQ6.6)/a NLO (MRSW) - 1. The 
dashed curves are the percentage variations from the central prediction between the upper and 
lower predictions obtained using the CTEQ6.6 PDF error sets. 



where gq is the Born cross section of Eq. devaluated with the 1— loop renormalization group 
improved value for and includes the full m& mass dependence of Eq. HI (Including the 
mj dependence has almost no numerical effect). 

The NLO QCD corrections are parameterized by the factor Aqcd and unlo is evaluated 
with the 2— loop renormalization group improved value for ~g~b(fi) 6 . The NLO QCD correc- 
tions to the bg — > bH process have been previously found in the S-ACOT scheme, which 



includes all effects of the finite b mass to 0(a s ). In the S-ACOT scheme 49|, effects of a 
non-zero b quark mass in the process bg — > bH are absorbed into the definition of the PDFs 
and to 0(a s ), we have schematically, 

&nlo,qcd(pp -> bH) mbjt0 = a LO (pp -> bH) mb=0 (l + 5 QC d) + 0(a 2 s ) 

olo^PP -> bH) mb=0 



*QCD 



°lo(pP -» bH) mb ^ 



Sqcd , (34) 



6 Our results agree with those obtained from MCFm[4^| and in Ref. [lj and are presented here only to 
facilitate comparison with A.wk ■ 
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FIG. 9: Tevatron results for the weak corrections to pp — > b(b)H with y/s = 1.96 GeV, > 
20 GeV, and | rjb |< 2. The solid black curve represents the contributions which cannot be 
factorized into an effective bbH vertex contribution and is always less than 1%. 

where 5 QC d is given in Ref. |fi|. Both the CTEQ and MRSW PDF sets employ the S-ACOT 
scheme and so our inclusion of b mass effects is consistent to 0(a s ). 

The QED and weak corrections are contained in Aqed and Awk, respectively. As 
discussed above, we do not present results for Aq ED , but assume they are negligible. The 
contribution of Awk results from the interference of the tree level amplitude with the 1-loop 
amplitudes shown above, which are generated numerically. We compare our exact results of 
Eq. |33]with an "Improved Born Approximation", IBA, which is obtained by replacing the 
tree level bbH vertex of Eq. [I] with the on-shell one loop electrpweak corrected vertex which 
can be found from the corrections to the decay H — >■ bb 34, 



V(H bb) = F(H bb) \l + Aqq D + Agf D + A% H K 
We define the Improved Born Approximation in an obvious fashion as 

a(bg bH) IBA (fi) = a(bg bH)M (l + A b ™ D + A% + A b ^j 



(35) 



(36) 



The IBA approximation assumes that the bulk of the weak corrections modify the bbH 
vertex. In the case of the SUSY QCD corrections to bg — > bH from squarks and gluinos, the 
Improved Born Approximation is an excellent approximation to the full rate [3]. 
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Results for the Tevatron are shown in Figs. [7] [8] and [9] The Tevatron plots have 
yfs = 1.96 TeV, | r)b |< 2.0 and require > 20 GeV. The NLO QCD corrections combine 
partons if AR = a/(A0) 2 + (At/) 2 < 0.4. For M H = 160 GeV, the scale uncertainty at NLO 
with a variation from fx = M H /2 to 2M H is ~ 10%, while for M H = 120 GeV it is ~ 8%. 
The PDF uncertainties are estimated in Fig. [8] where we compare the CTEQ6.6 predictions 
with those obtained using the MRSW2008 NLO PDFs 50], (with fi = M H /2), and find 
agreement between the 2 PDF sets to within better than 5%. The PDF uncertainties using 
the CTEQ6 error sets are also shown in Fig. [8] and are quite large, varying between 15% 
and 20% for the masses considered here 7 . 

Fig. [9] shows the size of the weak corrections as defined by Eq. [33] We note that the \i 
dependence of Awk is extremely small. The weak corrections are well approximated by the 
IB A of Eq. [36] (the dot- dashed line of Fig. [9]), with the remaining corrections (the solid line 
in Fig. [9] ) always less than 1%. Except near the W + W and ZZ resonances, A WK in the 
Standard Model is significantly smaller than the uncertainties from the QCD scale variation 
and the PDF uncertainties. 

At the LHC, we consider y/s = 7 TeV and y/s = 10 TeV, with \rj b \< 2.5, p b T > 25 GeV 
and AR > 0.4. The NLO QCD corrected cross sections are shown in Figs. [10] and [13] The 
NLO cross section is reduced by a factor of ~ 2.2 for M H = 150 GeV (with /i = M H /2) 
when going from yfe = 10 TeV to 7 TeV. At ^/s = 7 TeV, and M H = 150 GeV, the 
scale uncertainty at NLO with a variation from fi = Mh/2 to 2Mh is ~ 5%, while for 
Mh = 300 GeV it is ~ 9%. The PDF uncertainties for a/s = 10 TeV are estimated in Fig. 
H~4l where we compare the CTEQ6.6 predictions with those obtained using the MRSW2008 
NLO PDFs 50], (with /i = Mh/2), and find agreement between the 2 PDF sets to within 



better than 3%. The PDF uncertainties using the CTEQ6 error sets are also shown in 
Fig. [H] and are smaller than at the Tevatron, varying between 4% and 6% for the masses 
considered here. The PDF uncertainties are similar for ^/i = 7 TeV. 

The weak corrections are shown in Figs. [H] and [15] for Mh < 500 GeV. The IBA (Eq. 
36|) encapsulates the total weak corrections to better than 1% for Mh < 500 GeV. We show 
the weak effects for Mh > 500 GeV, along with the large Mh limit of Eq. [31] in Figs. [12 
and [16] For Mh = 1 TeV, the IBA underestimates the total weak corrections by about 



7 The PDF uncertainties obtained from the 40 CTEQ PDF error sets were previously obtained in Ref. [ij . 
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FIG. 10: Lowest order and NLO QCD results for pp -> b(b)HX at the LHC with y/s = 7 TeV, 
p\ > 25 GeV, | rjb |< 2.5, and Ar > .4. The renormalization/factorization scales are set equal to 
f i. 

3% at y/s = 10 TeV. For large Mh (Mh > 2M t ), the weak corrections are significant and 
are greater than 18% for Mh ~ 1 TeV. We note that the large Mh limit underestimates 
the weak corrections by about 5% at Mh = 1 TeV, implying that the log(M#) terms are 
numerically important. For heavy Higgs bosons, Mh > 500 GeV, the weak corrections are 
larger than uncertainties from PDFs and the scale choice, and it is meanful to include them 
in precision calculations. 

B. The m[, = Limit 



It is interesting to consider the — > limit of the bg — > bH amplitude. In this limit, the 
6-Higgs Yukawa coupling vanishes, = m^/v — > 0, and the tree level amplitude shown in 
Fig. [T]is identically zero. The first non-zero contributions to bg — > bH with m\, = arise from 
the squares of a subset of the 1— loop amplitudes shown in Figs. [5]and|6]and are 0(a s G 3 F ). 
The contributions which are non-zero in the m& — > limit involve the coupling of the Higgs 
to either a top quark or a pair of gauge bosons (and the corresponding Goldstone bosons). 



These contributions have been calculated in Ref. 



23| and we have checked that the squares 
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FIG. 11: LHC results for the weak corrections to pp — > b(b)H with yfs = 7 TeV, p\ > 25 GeV, 
and | r/f, |< 2.5. The solid black curve represents the contributions which cannot be factorized into 
an effective bbH vertex contribution and is less than 1% for Mjj < 500 GeV. 



of our 1— loop amplitudes reproduce their results in the = limit. Since these diagrams 
are not suppressed by a small b quark Yukawa coupling, they give a comparatively large 
contribution. At a/s = 7 TeV and Mn = 120 GeV, we find that the 0(a s G z F ) contribution 
with mj = is around 8% of the Born cross section shown in Fig{10] with our cuts. 

Although our calculations are purely Standard Model, we are, however, motivated by a 
very different scenario than the authors of Ref. 23|]. In models with an enhanced coupling 
of the b quark to a Higgs boson, the tree level amplitude can be significantly larger than in 
the Standard Model. In such models, it is important to understand the numerical effect of 
the interference of the tree level amplitude with the one-loop weak corrections. Future work 
will explore the role of the electroweak corrections in models with non-standard b quark 
Higgs Yukawa couplings, in particular the MSSM with large tan/3. 



IV. CONCLUSION 



We have computed the Standard Model weak corrections to the processes pp — > b(b)H at 
the LHC and pp — > b(b)H at the Tevatron. In both cases, the results are well approximated 
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FIG. 12: LHC results for the weak corrections to pp — > b(b)H with \fs = 7 TeV, p\ > 25 GeV, 
and | rjb |< 2.5. The solid black curve represents the contributions which cannot be factorized into 
an effective bbH vertex contribution. The dotted line is the large Higgs mass limit of Eq. [3TJ 



by including only the on-shell bbH vertex corrections, with the remaining weak corrections 
less than 1 — 2% for Mh < 500 GeV . This observation makes it straightforward to estimate 
the weak effects of non-Standard Model b quark Yukawa couplings on the bH production 
process. 

At the Tevatron, the weak effects are always much smaller than scale and PDF uncer- 
tainties and so can be neglected in the Standard Model. At the LHC, for large Mh the weak 
corrections can become significant and can be larger than scale and PDF uncertainties. At 
the LHC with y/s = 10 TeV, the corrections of O^y^J are ~ 18% for M H = 1 TeV. 
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